Probability Distribution of the Quality Factor of a 
Mode-Stirred Reverberation Chamber 
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Abstract — We derive a probability distribution, confidence 
intervals and statistics of the quality (Q) factor of an arbitrarily 
shaped mode-stirred reverberation chamber, based on ensemble 
distributions of the idealized random cavity field with assumed 
perfect stir efficiency. It is shown that Q exhibits a Fisher- 
Snedecor F-distribution whose degrees of freedom are governed 
by the number of simultaneously excited cavity modes per stir 
state. The most probable value of Q is between a fraction 2/9 
and 1 of its mean value, and between a fraction 4/9 and 1 of its 
asymptotic (composite Q) value. The arithmetic mean value is 
found to always exceed the values of all other theoretical metrics 
for centrality of Q. For a rectangular cavity, we retrieve the 
known asymptotic Q in the limit of highly overmoded regime. 



I. Introduction 

To date, the study of mode-tuned or mode-stirred rever- 
beration chambers (MT/MSRCs) - i.e., multi-mode cavity 
resonators furnished with a 'stirring' mechanism that produces 
time-varying excitation and/or boundary conditions - has 
mainly focused on the properties of the random electromag- 
netic (EM) field. Probability density functions (PDFs) for 
idealized and imperfect fields, including EM boundary-value 
problems [1], [2], were calculated and compared with mea- 
surements or simulations. A natural extension is the stochastic 
characterization of intrinsic EM parameters of instrumentation 
and devices subjected to random fields, e.g., wave and input 
impedances [3]-[6], antenna parameters [7], [8], etc. 

One of the fundamental parameters of a MT/MSRC is its 
quality (Q-) factor [8]— [21] . In the simplest model, Q is 
defined by a constant single value, as the ratio of the stir- 
averageqj stored energy (U) to the averaged dissipated power 

' The fact that the total energies stored and dissipated vary with changing stir 
state has been validated by experiments, which show that the measurement of 
Sn at a fixed frequency and source power P s exhibits random fluctuations 
with changing stir state. Therefore, the net forward power P^x = (1 — 
\Sn\ 2 )P B injected into the MT/MSRC fluctuates accordingly. 



(Pd), multiplied by the excitation frequenc^ uj [10]— [13] 
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The definition of this so-called 'effective' or 'composite' 
quality factor is inspired by the corresponding notion of modal 
Q for a single eigenmode of a static resonant cavity, i.e., 
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u, 



in np(j^mnp) 
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in which uj mnp , U mnp 



Pd,mnp i^mnp) 

and Pd,mnp take constant values for 
a selected mode specified by modal indices to, n, p. In a 
MT/MSRC, however, U and Pd fluctuate quasi-randomly as a 
function of stir state r. Hence, defining an instantaneous value 
of Q at each r as 



Q(u,t) = u(t)- 



(3) 



this Q is now a randomly fluctuating quantity with an associ- 
ated PDF /q(<z), correlation functions, etc., when considered 
across all r. For simplicity of notation, we shall further omit 
indicating the dependencies on u and t in ([3]). 

Compared to (fTJ, the definition ® is closer in spirit to 
the original concept of Q mnp - Firstly, OJ involves a ratio 
of quantities that exist physically at each t, as opposed to 
the formal ratio of mean values in (T|) that exist only in a 
r-averaged, i.e., mathematical sense. Secondly, because of 
propagation of uncertainties, any disregard for the random 
fluctuations of U and results in an underestimate of the 
level of fluctuation of other stochastic EM quantities that 
depend explicitly or implicitly on Q, in particular the standard 
deviations of the EM fields <J E ,{,) and cr„/(/) (cf. (l27b ), which 
are of fundamental importance. 

2 In general, the spectral power density g(u) and, hence, the spectrally 
averaged angular centre frequency ujq = J uig(w)dw/ J g(aj)daj vary 
as a function of stir state [22]. Therefore, spectral and ensemble averagings 
are strictly needed to replace w in (TJ by (uio). For narrowband excitation 
or nondispersive the fluctuations of ujo are usually negligibly small, 

whence ensemble and spectral averaging of ui can then be omitted. 
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An alternative but more restrictive approach to quantifying 
the uncertainty of Q was developed previously in [20]. There, 
second-order statistical characterization of Q was performed 
based on spectral and ensemble averaging of Q mnp for TE and 
TM eigenmodes. This permitted a calculation of the mean /iq 
and standard deviation uq for a rectangular cavity in which 
wall stirring causes modal perturbations. 

In the present paper, the use of sampled instantaneous (as 
opposed to averaged) values of Q allows for the calculation of 
the complete PDF /q(<z). This provides a more comprehensive 
characterization compared to mere first- and second-order 
moments. For simplicity, the analysis is based on ensemble 
distributions of U and Pd for ideal Gaussian EM fields, 
as opposed to their sampling distributions [23], [24]. This 
implicitly assumes that a sufficiently large and theoretically 
infinite number (N) of statistically independent stir states for 
the field is generated by the stir process (N —> +oo). It 
will be found that fq(q) then satisfies a Fisher-Snedecor F- 
distribution, whose two numbers of degrees of freedom (DoF) 
both depend on the number (M) of simultaneously excited 
cavity modes per stir state. When M increases, fq(q) evolves 
from a positively skewed PDF for low M toward a Gaussian 
(normal) PDF, accompanied by a reduction in mean value and 
absolute or relative spread of Q. 

The results apply generally to cavities with arbitrary ge- 
ometries, including irregular shapes, but will be illustrated 
with explicit expressions for simple (integrable) rectangular 
cavities. Except in Sec. IIII-DI we do not distinguish between 
mode-tuned and mode-stirred methods of operation, insofar 
as only quasi-static fields are considered. An exp(jwt) time 
dependence is assumed and suppressed throughout. Different 
types of averaging will be performed: we shall use the 
notations (-)v, (-)s and (•) to represent spatial averaging 
with respect to the volume V, surface area S, and ensemble 
averaging with respect to cavity stir states r, respectively. The 
ensemble average assumes equal V and S throughout, as a 
prerequisite for constant average spectral mode density, and 
only involves perturbations of shape or aspect ratio(s) across 
different realizations. 



II. Unstirred chambers with single-mode 
excitation: deterministic U and P d 

For the purpose of establishing definitions and extending 
the analysis later to stirred multi-mode cavities, we briefly 
review some basic results for the Q of a single eigenmode 
in a static (unstirred) cavity at its fixed angular resonance 
frequency u! mnp . In this case, the modal Q m n P is single- 
valued. In an unstirred cavity, the local electric and magnetic 
modal amplitudes Eq and Hq at any location r inside V are 
time invariant. For a lossless linear time-invariant isotropic 
homogeneous medium filling the cavity interior, D_ = €qE_ and 
P = /j,qH,, whence the electric and magnetic stored energies 
U c and U m are proportional to the spatial integrals of the local 
intensities |P(V)| 2 and \H_(r)\ 2 , respectively. For steady state 
excitation, U e = U m to first approximation at sufficiently high 
frequencies. The total stored energy U = U c + U m can thus 



be expressed as 
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where an asterisk denotes complex conjugation. 

To obtain a corresponding expression for Pd, the relevant 
quantity is the tangential magnetic field P t at a location 
r s on the cavity's interior boundary surface S — dV with 
unit local inward surface normal l n (r s ). Conduction loss 
in the wall yields a nonvanishing tangential electric field 
E_ t = P w (P t x lj ^ at S that can be envisaged as 
a surface layer of magnetic current produced by H_ t as an 
equivalent boundary source and dissipated by S, in addition 
to the sheet of surface charge produced by the normal electric 
field in a lossless cavity. The time-averaged absorbed energy 
is the spatially integrated real part of the normal component 
of the local Poynting vector, S n = Re[l„ • {E x H*)]/2, i.e., 



Pd = 



^ / \UL s )xH(Ls)\ 2 dS 
J s 

[ \H t (r s )\ 2 dS 
Js 



2a w S 



2cr w 5\ 
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in which P w = 1/(ct w <5 w ) = ^/w^ w / (2cr w ) represents the 
per-unit area surface resistance of the interior cavity wall, 
S w = y/2/ (w^ w (T w ) is its skin depth, and where /x w = n WyY ^o 
and er w are its permeability and conductivity, respectively. 
Substituting © and © into © yields [25, sec. 10.4] 



Q 
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If only one resonant mode is excited, then the local field 
amplitudes throughout V are characterized by a single modal 
amplitude value P m „ Pi o, whence © can then be written as 



Q 



IP 



mnp,0 | 



Iv ^ it) 



2 dV 



>^\ 2 fs\t mnp (r_ s WdS> 



(7) 



where 



^ is the real-valued magnetic eigenvector, and 
Hmnp,{t,)o i s i ts associated complex-valued amplitude. 

As an example, consider a rectangular cavity with V = 
l x l y i z with a single excited mode whose local amplitude 
\H_ m (r)\ at r = xl x + yl + zl z can be expressed as 



IP 



in up 



(r)l 



|^ | J sin(/L mn p ja; x) 1 
\ sin(fc mn p } y2/) J cos{k rnn p^ z z) 



A similar expression for the electric field E_ mnp {r) with ampli- 
tude \E mnP} o\ applies. Thus, the amplitudes are unmodulated 
(sinusoidal or constant) with respect to r, in all directions. 
For a transverse mode, any valid combination in ([8]) consists 
of two spatial harmonic functions along two orthogonal di- 
rections l a and lp combined with the unit function in the 
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third direction 1 (i.e., cos(fc mnPj7 7) = 1), where a, /3, 7 £ 
{x, y, z}. However, in MT/MSRCs, we are primarily interested 
in overmoded conditions at high frequencies, where the vast 
majority of modes are hybrid. For a hybrid mode satisfying 
the EM boundary conditions, a valid combination in (0 is 
the product of three such harmonic functions. Substituting 
(0 into (0 and (0, together with J " sm 2 (k mnp ^ a a)da = 
J a cos 2 (k mnpa a)da = £ a /2, we obtain for a rectangular 
cavity with conducting boundaries 



16 



mnp,0 I 



5" / ^W , |2 

Pi — — \j \-tlmnp,t,0\ , 



(9) 



in which the difference between the factors 1/16 and 1/8 
results from the fact that a hybrid mode generally exhibits 
three magnetic field components for the interior field, whereas 
only two nonvanishing tangential components of this field exist 
on the surface. For a transverse mode, the factors 1/16 and 
1/8 in © are replaced by 1/8 and 1/4 (or possibly 1/2, but 
with vanishingly small contribution, when the cavity surface is 
locally perpendicular to the transverse direction of the mode), 
respectively. From (0 and (0, it follows that 



Q = 



v 



mnp,0 I 



(10) 



For any single mode in a rectangular unstirred cavity, 
Hmnpfl and P„ mPi t,o are constant with respect to location 
and time. In nonrectangular (e.g., cylindrical) cavities, the 
eigenmodes no longer consist of spatial harmonics. For non- 
separable geometries, they may not even exist in closed-form 
expressions. Hence the ratio (\H^(r)\ 2 )v / (\H_ t (r s )\ 2 ) s is in 
general different from \H mnpfi (r)\ 2 / {2\H mnp . t fi(r s )\ 2 ), in 
which case < TT~0b is then multiplied by a shape dependent factor 
h (cf. (El). 

III. Stirred chambers: quasi-random U and P d 
A. Single vs. multimode excitation 

When invoking mode tuning or mode stirring, the eigen- 
modes and therefore Q evolve with stir state r. The charac- 
terization in Sec. [TT] for a single mode is then only meaningful 
in a statistical sense. We denote {P(i"|r)} and {H.( t )( T \L(s))} 
to represent ensembles of sample sets of N stir states of the 
field at arbitrary r_( s) ■ These stirred local fields vary randomly 
in spatial orientation, magnitude and phase as a function 
of t. Generally, P(r|r) and H_ t (r\r s ) have three and two 
nonzero complex-valued (in-phase and quadrature) rectangular 
components H a (a = x, y, z), respectively. Hence, if at any r 
only one cavity mode is excited that is randomly perturbed by 
the stir process (i.e., random single-mode excitation), then the 
local fields H_ (r|r) and H. t (r\r_ s ) in (0-(0 are now random 
processes with six and four DoF, respectively. 

In practical (non-superconducting) overmoded MT/MSRCs, 
significant multimode excitation occurs typically. Spectral 
overlap of nondegenerate modes having nonzero absorption 
bandwidths causes intermodal coupling, such that even a 
single-frequency (CW) source then excites simultaneously 
several modes with different Lu mnp . Alternatively, the spec- 
trum of a wide-band source may encompass several u) mnp 



of the MT/MSRC, so that corresponding (non)overlapping 
modes may be simultaneously excited (e.g., multitones in a 
multimode laser or in certain communications protocols). Let 
the number of simultaneously excited modes per stir state 
be denoted by M . A physical estimation of the value of M 
is given in Appendix ??. For arbitrary r, the cavity field is 
the resultant of the weighted superposition of M participating 
modes (random walk model), producing a spatial modulation 
of the amplitudes Eo and H(t.)o across the cavity's interior. If 
the structure of the cavity is sufficiently irregular at cj, then 
this spatial distribution is quasi-random. This (static) spatial 
variation is additional to the (dynamic) fluctuations of the 
local field caused by stirring and yields a 4-D spatio-temporal 
random field. For spatially random fields, (0 and (0 remain 
valid in a statistical sense, i.e., at arbitrary r. Thus, for a wall- 
stirred rectangular cavity, we now have instead of (l9l-(fT0li. 



U(r) 
Q(t) 



16 

S /w/i 



Mr) 



2(7, 



-(l^o(r)| 5 



V (\H Q (r)\ 2 )v 



(11) 
(12) 

(13) 



Mwj-(5wS (\H t:0 (T)\ 2 )s' 
Each individual mode acts as a 'channel' for storage and 
dissipation of energy and increases by one unit the number of 
ways in which the value of each Cartesian component of the 
resultant field phasor can be obtained. The DoFs of the spa- 
tially integrated field intensities leading to U and Pd increase 
accordingly. For each mode, E_ mnp and H_ mnp are physically 
(i.e., deterministically) related via a wave impedance dyadic. 
Following (0-01, the increase of the number of DoF of U by 
each mode is therefore the same as for U c and U m individually, 
i.e., six, whereas the corresponding increase for Pd is four. 

Finally, assuming that the stirring process is sufficiently 
efficient to be capable of generating a very large (theoretically 
infinite) value of N that produce independent and identically 
distributed ideal Gaussian Hu) (i- e -' N — > +00) and as- 
suming that wide-sense ergodicity of the fields holds (i.e., 

VHM-r) = Vh{t\t) and °i£(r|T) = °ir(T|r))> such that the 
spatial distributions of each H^ t ) a is identical to its ensemble 
(i.e., stir) distribution, U then exhibits approximate^ a x^m 
PDF, and Pd has approximately a x\m PDF across the stir 
states and cavity interior, i.e., (l34l and (f35T > hold with r — 3M 
and s = 2M, respectively. 

In practice, the values of M and N evaluated for different 
stir processes are often strongly positively correlated. For small 
M, the practical stirring performance may then be significantly 
compromized. If the assumption N —> +00 becomes unsus- 
tainable, then the use of Bessel K sampling distributions for 

3 The incoherent superposition of M participating modes with equal x| 
energy distributions presumes that energy is equally partitioned across these 
modes. For overlapping modes, this is only approximately true because the 
partitioning depends on the source's spatial location and on the specific 
spectral distances of the u) m i n i p i relative to the excitation frequency uj. They 
should be Lorentz weighted accordingly in the superposition. 

4 On S, the nonorthogonality of irrotational and solenoidal magnetic eigen- 
vectors (i.e., solutions associated with boundary conditions of magnetic type) 
causes the total power loss to deviate from the sum of the power losses of 
individual modes [25, Sec. 10.4]. This nonorthogonality results in the DoF 
2s to be somewhat less than AM. 
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S.[t)> U and Pd\U [24] offers an appropriate framework for 
characterizing /q(<z) [27]. In practice, the case M = 1 often 
(although not exclusively) involves excitation at wavelengths 
that are not small compared to cavity dimensions, while also 
the modal overlap is small. In this case, the stir process is 
also typically (but not inevitably) less efficient, whence the 
Xq and x\ PDFs of U(r\r_) and Pd(T~\r s )\U (r|r) and, hence, 
the PDF < n~4T > are then only approximately valid. 

Although we shall further focus on the case where M 
is a constant integer with respect to stir state, one may 
envisage a situation where its value could fluctuate as a 
function of r. In this case, we can estimate M by its mean 
value T,t=i M ( t )/ N > which 

may be fractional, yielding a 
generalization of x\m an d PDFs to gamma PDFs. Values 
of M smaller than unity represent the case where, on average, 
less than one mode per stir state is being excited. 



B. Probability density function and statistics of Q 

Regarding the relationship between U and Pj, the boundary 
field H t (r s ) is deterministically related to the interior H(r), 
because of field continuity and EM boundary conditions. 
Therefore, U and Pd are not statistically independent. Never- 
theless, their joint PDF can always be expressed as the product 
of the marginal PDF of U and the conditional PDF of P<j 
given U, viz., fu,p d (u,p d ) = fu{u)fp d \u(u,Pd\u). Based on 
this factorization, /q(<?) is derived in Appendix ?? as the ratio 
of a x\r distributed ljU and a x\ s distributed Pd, for general 
values of r and s, resulting in a Fisher-Snedecor F-distribution 
with (2r, 2s) DoF (cf. eqs. ((37), gO), dHJ, © and (@7)). 
Assigning r = 3M and s = 2M, the PDF of Q is then 
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B(3M,2M) ( 9 . 
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(14) 



valid for M > 1/2, representing a F(6M, 4M) PDF. Figure Q] 
shows ( fT4b for selected values of M. For M ~ 1, significant 
positive skewness and smaller kurtosis are observed, whereas 
for M — > +oo, the PDF evolves to Gaussian normality. 
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Fig. 1. Scaled PDF /q(q) of normalized Q (i.e., in units (Q) — 
[3M/(2M - l)][(/iVy(/i w , r 5 w S)]) for selected values of M. 

The (arithmetic) mean value, standard deviation and coeffi- 



cient of variation for ( TBI follow from (|43T>— (f45l> as 

3M h V 
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(M > 1) (16) 
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- 1, (M > 1). 



3M(2M - 2) 
Their corresponding limit expressions for M 3> 1 are 

1 \ hV 



3" 



2M J fi w>T S w S 
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(17) 

(18) 
(19) 
(20) 



The dependencies of (fT~5l> — (f2Qb on M are shown in Fig. [2] For 
M — > +oo, the mean (Q) reduces asymptotically to half its 
value for M = 1, i.e., to 

* 3 ^ (21) 



The residual mean A(Q)/Q oc = ((Q) - Qoo)/Qoo = 
1/(2M — 1) is positive and asymptotically inversely pro- 
portional to M. For h = 1, the result (fJTJ was previously 
obtained [10], [13], [20], [26], whereas the finding that (Q) = 
2Qoo when M — 1 agrees with the findings in [15]. This 
demonstrates that consistent asymptotic results are retrieved. 
In Sec. IIII-CI (Q) will be compared to other measures of 
centrality for Q. For M = 1, the <jq and vq are undefined 
(cQ; V Q ~^ whereas for M — > +oo they asymptotically 

approach zero proportionally to i.e., more slowly than 

A(Q) /Qoo. In summary, for M — > +oo, the limit PDF of Q 
is a normal distribution M{hq, ctq), i.e., 




2m) 0o °' V 6M 



(22) 



Parenthetically, if M ^> 1 then the ratio of the averages 
(ojU) and (Pd\U) is substantially different from the averaged 
ratio (uU/(P d \U)). Indeed, for M = 1 with ([Toll, (B]> and 
(EH), i.e., (|P | 2 ) = 6<r^, (0 for a x| distributed |P | 2 while 
(IH^o) 2 ) = ^ a \i(i) f° r a xl distributed |P t . | 2 , we arrive at 



>U) 



hV 



3hV 



(Pd\U) 



1 / ujU 
2\PJU 



l 2 > 
(M 



2/z, 
= 1). 



(23) 
(24) 



In view of d2Tb . this result shows that replacing (toU/ (Pd\U)} 
by (luU) / (Pd\U) as in (fl} is an approximation, but justifiable 
when M 3> 1, e.g., in overmoded regime. In fact, comparing 
(03, <ED and d23l using (l39l l for general M shows that Q c s = 
Qoo for any M. Thus, the definition of composite Q c ^ neglects 
the effect of M on the fluctuation and value of Q. 



5 



10 



10' 



10 u 



10" 



10" 



10 



10 



M 



I I 

■ v. ft 


• A<Q)/Q 
V Q 





10 



"D 

0) 

E 



10 



10 



O 10" 



E 



10" 




20 40 60 

(a) 



80 100 



Fig. 2. Mean, residual mean, standard deviation and coefficient of 
variation of Q as a function of M, normalized by Qoo- 



As an alternative to <tq, the spread of Q can also be 
expressed by an 77%-confidence interval for Q. For a chosen 
confidence level 77, the boundaries qg and q u of this interval 
are calculated by inverting the cumulative distribution function 
(CDF) d48l>— (f50b. i.e., by numerically solving 

F Q {q Lu ) = 1 - l u>u (2M, 3M) = 1±(?7 2 /100) , (25) 
where 7^ u (•, •) is a regularized incomplete beta function with 



3M 



ft, 



2M - 1 (Q) 



(26) 



Figures [3^ and show these boundaries normalized by the 
median of Q (cf. (f30b— (f3Tb ~) or by (Q), as a function of 77 or 
M, respectively. The interval width rapidly increases with 77 
most prominently when M ~ 1 and 77 > 90. For larger M, 
the spread is considerably reduced because of the effect of 
aggregation of modes, which can also be achieved through 
multiple stirring mechanisms, multiple sources (antennas), 
increased EM losses, etc. For M = 1, 3, 10 and 100, 
the 95%-confidence intervals for Q/(Q) are [0.080,4.599], 
[0.301,2.590], [0.545,1.713] and [0.833, 1.192], respectively. 

The standard deviations of the complex-valued analytic EM 
fields E = E'- ]E?' and H = H' - jH" can be estimated 
on the premise that, in steady state, the dissipated power 
equals the transmitted power, i.e., P<i(t) = Ptx(t), due to 
conservation of energy. With <(3j and d39l for n = s = 2M, it 
follows that (P d \U) = MSa 2 H , {l) /2, whence 
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(27) 



where 770 = yj Jj® Jeo is the stir averagecH input impedance 
of the MT/MSRC. Thus, like for Q, the standard deviation 
of the stirred EM field decreases proportionally to 1/yM. 
This is a result of intrinsic averaging of fields caused by the 



More accurate estimates for a 



T 2 

rating the random fluctuations of the input impedance dyadic Z_(t) [3]-[6] 
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are obtained by incorpo- 
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Fig. 3. Upper (q u ) and lower (qi) boundaries of 77%-confidence 
intervals of Q: (a) normalized by med(Q), as a function of confidence 
level rj%, for selected values of M; (b) normalized by (Q), as a 
function of M, for 95%, 99% and 99.5% confidence levels. 



simultaneous excitation of modes. For the total (3-D vector) 



fields, cr|, w 



3cr 



and a 2 H , (l) 



3ct 2 , M . 
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C. Other measures of location for Q 

Because of the primary practical interest in the central value 
of Q, we explore a few other measures of location (centrality) 
as alternatives to the arithmetic mean (Q). Compared to such 
other metrics, (Q) represents the 'centre of mass' of the PDF 
and minimizes the expected mean squared deviation of the 
sample values of Q. The (Q) is known to provide the most 
stable measure of centrality when comparing values obtained 
from different sample sets of data. However, it is not the 
optimal measure of centrality for ensemble data, particularly 
when the PDF is significantly skewed, as in the case of 
relatively small M. In the latter case, the mode (for unimodal 
data) and the median are more representative parameters. 

The generalized mean (Q) a = (J^° q a fg(q)dq) 1 ^ a can be 
calculated with the aid of d42l and (l2Tb as 



_ 2 / r(3M + ffl)r(2iir-a) 

{Q)a - 3 y r(3M)r(2M) 



1/0 



Qc 



(28) 
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where a is a chosen real parameter. The particular cases 
a = — 1, 0, 1 and 2 correspond to the harmonic, geometric, 
arithmetic and RMS averages, respectively. Figure |4] shows 
(Q) a as a function of a for selected values of M, after 
normalization with respect to (Q) = (Q)\. For any a < 1, 
(Q)a is always smaller than (Q), a fortiori for small M. For 
M — > +oo, all (Q) a merge to (Q) irrespective of a. 
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mcd(Q)/(Q) increases from over 0.5308 to 1, whereas 
mcd(Q)/Qoo decreases from 1.135 over 1.062 to 1. 

In practice, experimentally determined values of Q are 
nearly always reported to be considerably smaller than Q c ff = 
Qoo, typically by a factor 0.2 or 0.5 to 1 [7], [9], [28]. On 
account of (l32l . choosing mod(Q) or a fortiori (Q)c a <-2) as 
an a priori theoretical estimate may offer better ad hoc quan- 
titative agreement than (Q) and mcd(Q). Physically, however, 
additional loss mechanisms are at the root of lowering Q [10]. 
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Fig. 4. Generalized mean (Q) a normalized by arithmetic mean Fig. 5. Mode-to-mean and median-to-mean ratios mod(Q)/{Q), 
(Q) = (Q)i as a function of a, for selected values of M. med(Q)/(Q), mod(Q)/Qoo and med(Q)/Q ao as functions of M. 



The statistical mode is the most probable (or most frequent) 
value among the values of the population (or sample data set), 
i.e., for {Q(t)} across one rotation of a mode stirrer. Solving 
d[fQ(q)]/dq = using ([T4l >. the mode-to-mean ratio is 



mod(Q) 



1 



1 

3M 



1 



1 



3M- (29) 



(Q) V~ 3M J V 2M + 1 / 
Another measure of centrality is the median, defined as 
med(Q) = Fq 1 (0.5) and obtained from d48l-d50l by solving 

1 

2 : 



I u (2M, 3M) 



(30) 



where 



1 



3M 



2M- 1 



ied(Q) V 

(Q) J 



(31) 



A numerical approximation of med(Q) is obtained from 
d28l for a — —0.065. Unlike (Q), the median minimizes 
the expected absolute deviation. It is a robust measure of 
centrality, being less sensitive to the shape of fq(q). This 
is particularly attractive because of the cited difficulties of 
characterizing the precise PDFs of U and Pd when M > 1. 
Comparing these metrics, for any M > 1/2, the ordering 



mod(Q) <Qoo< med(Q) < (Q) < Q RMS 



(32) 



applies, together with mod(Q) < (Q) a < Qoo when — 1 < 
a < -0.185, Qoc < (Q)a < med(<3) when -0.185 < a < 
-0.065, and med(Q) < (Q) a < (Q) when -0.065 < a < 1. 

Figure [5] shows that the ratio mod(Q)/(Q) increases from 
when M — s- 1/2, over 2/9 at M = 1, to 1 when M +oo. 
With the same marker values of M, this corresponds to 
mod(Q)/<5 00 increasing from 1/6 over 4/9 to 1. The ratio 



D. Q in mode-stirred vs. mode-tuned chambers 

The previous analysis assumed quasi-static operation, 
through sufficiently slow stepping or scanning (mechanical or 
electronic). General considerations of mode-stirred vs. mode- 
tuned operation focusing on acquisition time and nonstationary 
effects were given in [21], [22], [29]-[32]. Here, we comment 
on aspects of mode stirring relating to Q and its PDF, and we 
restrict ourselves merely to some general remarks. 

In quasi-stationary MT/MSRCs, the rate of change of the 
cavity field between stir states is small compared to the 
rate of energy fill and dissipation, whence the values of 
U(r,r) and Pd(r_si T ) remain unaffected by this rate. Quasi- 
stationary mode stirring may results in purely local uniform 
temporal averaging of the fields across an interval of stir 
states [0,7] varying with time. This averaging does obviously 
not affect the spatial averagings of U(r\r) and Pd(rsl r )' 
yielding again dTTb and ( TT2b . Ensemble averaging of U (r) and 
P<j(t) across [0,7] does not affect (U(t)} and {P a {t)\U (t)) , 
whereas both cr^,% and cr^,,^,-. scale by the same variance 
function 77- in quasi-stationary conditions [22]. Consequently, 
the ratio au/ap A \u and hence /q((j) remain unaffected by 
local averaging. However, the existing coupling between TE 
and TM modes caused by ohmic losses is further increased 
by continuous rotation, albeit as a second-order effect of 
velocity of rotation [33], [34]. In turn, this has a positive but 
second-order effect on Pd- Recent measurements of (Q) in 
an overmoded MSRC [35] appear to support these findings, 
including a marginal but systematic decrease of (Q) observed 
for increased stir speed using the mean power approach. 

For nonstationary MT/MSRCs, the situation is more in- 
tricate. Expanding wave fronts emanating from the source 
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impinge and get partially absorbed by the cavity walls, after 
traversing the cavity interior. Therefore, the cavity fill and 
dissipation processes per stir state are not synchronized and 
may be affected differently during nonstationary stirring in a 
resonant environment, as boundary and excitation conditions 
vary in a rapid manner. This situation is to be avoided, in view 
of the definition ([3} requiring matching (i.e., comparable) pairs 
of U(t) with P<j (t). Hence, for a nonstationary MSRC, one 
may need recourse to the physically less meaningful definition 
((T), thereby using nonuniform weighting [36]. 

IV. Conclusion 

In this paper, we derived a PDF for the Q of a MT/MSRC, 
with the number of simultaneously excited cavity modes (M) 
as a distribution parameter. The analysis assumed ideal \ 2 
ensemble distributions for stored energy and dissipated power, 
associated with unbiased circular Gaussian distribution of the 
stirred cavity field, in space and time (stir domain). Perfectly 
efficient mode stirring, i.e., a theoretically unlimited number 
of statistically independent stir states of the local field (N — » 
+oo), was also assumed but can be relaxed by using sampling 
distributions to replace the ensemble distributions. In practice, 
this condition on N is usually more closely achieved at short 
wavelengths (A -C V 1 ' 3 ). Also, strong positive correlation typ- 
ically exists between M on N in practical stirring techniques, 
which may leave the expressions for /q(q) to be somewhat 
approximate in practice when M ~ 1. 

With these idealizations, Q was found to exhibit a Fisher- 
Snedecor F(6M, 4M) PDF, given by (O. For general M, its 
arithmetic mean value (Q), standard deviation and coefficient 
of variation were obtained in (TT3T> — (TTTb. For M — > +oo, the 
PDF approaches a Gauss normal limit PDF (122) . Confidence 
intervals for Q were calculated numerically from (I25ll-(l26ll 
and indicate that the spread of Q can be considerable, even 
when M ^> 1. Expressions for the standard deviations of 
the EM fields were obtained in (l27l i. Alternative measures 
of centrality (viz., the generalized mean (1281 1. mode (|29l and 
median (f30l> — (l3TT>"> rank according to (I32t . producing values 
that are always smaller than (Q), for any M but especially for 
M ~ 1. Specifically, the most probable and most robust central 
values of Q are smaller than (Q) by factors ranging from 2/9 
and 0.5308 (for M = 1) to 1 (for M +oo), respectively. 
Compared to the mean and especially the median, the mode 
and generalized mean with index a < — 2 agree quantitatively 
with typical measured values of Q that are much lower than 
a priori estimated (Q). 

Regarding extensions of the present analysis, the nonorthog- 
onality of the magnetic eigenvectors on the boundary in case 
of multi-mode excitation is known to affect the total dissipated 
power. This requires characterization and quantification in 
the context of dynamic cavities. Generalization to nonzero 
excitation bandwidths is required to quantify the resulting 
increase of the uncertainty of Q. The effect of the non- 
integrability ('complexity') of real MT/MSRC enclosures - 
including the geometry of a mode stirrer - on h in d58l l 
and the effect of elongation or flattening of the cavity shape 
on /q(<?) deserve further attention. A rigorous generalization 



to fluctuating M during stirring (cf. Sec. IIII-AI) . caused by 
entrance and exit of modes in the cavity bandwidth [37], 
is of interest. The PDF /m(w) is expected to depend on 
the eigenfrequency spacing statistics [4] and their dynamics, 
which, in turn, depend on the integrability of the cavity shape. 
The additional uncertainty caused by fluctuating M is expected 
to increase the width of the confidence interval of Q. Finally, 
the (auto)correlation function of pq(At) provides insight into 
the rate of fluctuation dQ/dr, i.e., the stir sensitivity of Q(t). 
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Appendix 



A 



The PDF of Q = toll /Pd can be calculated from the joint 

PDlf] / u[ /,p d (wu,pd) of loU and P d = loU/Q, as 



1 f + °° 

/q(<?) = — / Uu, 
1 Jo 



LOU . 

;■ ( urn, — ) \wu\ d(cju), (33) 



in which f u u,p d (wu, pa) = fu,u(uu)fp d \ u u(j?d\uu), with 
= fu( u — £// w )/M f° r the auxiliary variate Y = loU 
and fp^uuiPdluu) = f Pd \u{Pd\u) for deterministic ui. For 
ideal Gaussian interior and surface fields, the associated energy 
density U and conditional dissipated power Pd | U exhibit xlr 
and xL PDFs, respectively, i.e., in self-sufficient form [22] 



r r/2 



s s/2 



r— 1 

u exp 



(34) 



-p d s 1 exp 



Pd ,05) 



fu{u) 

r ( s ) cr p d \u V a P d \u 

with which (f33j reduces, with the aid of [38, (3.351.3)], to 

/q(?) = ^ J fu(u)fp d \u [~J u du ^ 36 ) 
s ujajj 



B(r,s) VV r o-p d \u 



8 UJ<7U 



r+s ' 



(37) 



(38) 



where the complete beta function is calculated as 

Wr s) = T{s)T{r) = {S - l)] (r ' 1)! 
V ' ' T(s + r) (s + r-1)! ' 

Alternatively, (|37| | can be expressed in terms of average values 
serving as distribution parameters, because any X2n distributed 
variate X - being the sum of squares of 2n independent 
and identically distributed Af(0, a) variates - has its standard 
deviation and mean value related via [23] 



(X) _ r 2 

as a result of which d37b can be written as 



/«(?) 



s lj(U) 
B(r,a) \r(P d \U) 



r (P d \U) 



r+s ' 



(39) 



(40) 



6 It is assumed here that the excitation is sufficiently narrowband so that ui 
can be considered to be a deterministic (constant) quantity. 
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Equations d37T i and (|4Qb represent a Fisher-Snedecor F- 
distribution with (2r, 2s) degrees of freedom. 

Now we wish to express the ratio of the field-dependent 
distribution parameters crjj/crp d \u or (U) / (Pd\U) in d37| i or 
(l40l in terms of field-independent cavity design parameters. 
The evaluation in Appendix ?? yields these ratios as 
Hence, (f37l > and (t40b can be written as 



fab) 



i 



hV 



B(r,s) \fJ. w , T 5 w S 



hV \ 



r+s ' 



(41) 



in which ft, is the average shape factor of the stirred cavity, 
defined by 458}. 

With the moments of Q of order i defined as (Q l ) = 
Jo°° Q l fQ{q)dq and evaluated using [38, (3.194.3)] as 



T(r + i)T(s-i) ( hV 



(42) 



r(r)r(s) y^w.rSwS, 

the arithmetic mean, the standard deviation and the coefficient 
of variation are obtained as 

r hV 



/'Q 



(Q) 



S-l Mw,r^w<S'' 



(s > 1) 



(43) 



V(Q 2 ) W 



r(r + 1) 



hV 



(s — l)(s— 2) \s— 1/ ^ Wjr (5 w 5" 

(s>2) (44) 



f (s - l)(r + 1) 
(s - 2)r 



-1, (s>2) (45) 



respectively. The form fiTl i can therefore be rewritten in self- 
sufficient format in terms of (Q), i.e., 



f Q (q) 



i 



s - 1 



valid for s > 1, or alternatively in terms of erg, as 
1 



r-f-s 



(46) 



/«(«) 



( s -l)2( s _2) 



B(r, s) \ y r(r + l)(s - 1) - r 2 (s - 2) ' 
x q - ^ , (47) 



/ (s-l) 2 (s-2) 
9+ A/ r(r+l)(s-l)-r 2 (s-2) fJ Q 



r+s ' 



valid for s > 2. Both d46b and d4Tb depend neither on field 
parameters nor on cavity design parameters. 

The CDF of Q corresponding to the form d46*l ) is 



F Q (q) = l-I^s,r), (s>l) 



where 



A 1 



^- 1 (l-t) r - 1 dt 



B(s,r) J 

is the regularized incomplete beta function with 

r q 



-i(Q) 



(48) 



(49) 



(50) 



Appendix 

We consider statistically independent, circular Gauss nor- 
mal, analytical fields E and H_ . From ©, © and <[8j>, 

u(u) ^{i v m^dv) 

( Pd \U) l 2 ^ {Js \H t (r s WdS) 

j 2u, m a (J v \H(D\ 2 dV) 
V Mw, r (/ s |ff t (r5)l 2 d5)' 

For general multimode excitation, 



(51) 



(M) 

Hmnpfl 0, 



■nm p 



(r) 



(52) 



where X^mni denotes a sum over M simultaneously excited 



modes. Upon ensemble averaging, the cross product terms 
(incoherent field terms) in the summation vanish and we obtain 

(M) \ 



(Xiafc)iw) 




(X 









(\H \ 2 1 £iW^i 2dy ) (53) 

\ ^ mnp I 

M (\H \ 2 ) (J v m L )\ 2 dVy (54) 



The approximation in ( 1531 arises as a result of replacing the 
superposition of modal amplitudes H mnP: Q by a single random 
amplitude Hq, for arbitrary r. The equality d54"i l holds provided 
that the output impedance of the source is matched to the 
input impedance of the cavity for each stir state r, such that 
each source amplitude H mnp ,o is independent of r and, hence, 
independent of the variation of <b as a function of r. The 

—mnp 

modal triplets can be omitted when the integral of \d> I 2 is 

r ° '—mnp 1 

independent of k m np (cf. Sec. HJl. Similarly, 

(/j^fcs)| 2 d^ = M (\H t ^) (Jj£(r s )fdSy (55) 
Substituting d54l . (I5"5l and (O into d5TJ yields 

7 w([7) /f fey 



°>d|E/ 



r (P d \U) V s Vw.rdwS 
r 2 (/ y |<Kr)|W) 



where 



S /iw,r<$w (/ s |0(rs)| 2 dS') ' 

2 ((I0(£)| 2 )y) 
«l<Mlls)l 2 >S> 



(56) 
(57) 

(58) 



defines an average geometry (shape) factor of the stirred cavity 
that maintains constant V and S. For a wall-stirred rectangular 
cavity, |</>(x)| 2 = sin 2 (k x x) sin 2 (k y y) sin 2 (k z z) whence h 
1 and d56l) becomes 



UJCFlJ 



V 



a P d \U V 2 /J, w rS^S 



(59) 
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Appendix 

For practical application, an a priori estimate of the number 
of simultaneously excited cavity modes, M, based on cavity 
design parameters is of interest. 

The average mode density inside a cavity of volume V 
and surface S operated at CW frequency / is given by the 
generalized Weyl density [39], [40] as 

2 



d/ 



SttV (J 



c \ c 

1 

07TC 



4 

37TC 



d.s 



lav q{l) 
[tt — <^(n)][7r — 5<^(r)]dZ 



(60) 



ids ¥>(zi) 

up to terms of order (//c)~ 2 or smaller, where M c represents 
the cumulative spectral mode count from dc to /. The other 
symbols in ( f60b were defined in [40]. For A <C V 1 ^ 3 , i.e., 
V(f /c) 3 ^> 1, the number of modes SM C inside a narrow band 
of width Sf = (f/Q(f)) — f/(Q(f)) can be approximated 
by retaining only the first term in ( f60b and using d43l to yield 

SM c (f)^^^5f = 8n S -^^^(l-) 3 . (61) 
df r h \cj 

A CW source excite^ all overlapping modes within Sf 
simultaneously, hence SM C = M, With r = 3M and s = 2M, 
(f5TT > yields M as a solution of the quadratic 



where 



3M 2 - 26M + 6 = 0, 



s (f 

b — 8tt 



h 



(62) 



(63) 



On physical grounds, the solution with the negative sign is 
discarded because M is known to increase with /. Therefore, 



M(f) 



b + Vb 2 - 36 



whose asymptotic expression for / 

MooU) = — 



-oo is 26/3, i.e.. 



Qo 



(64) 



(65) 



i.e., proportional to / 5 / 2 and a^ 1 ^ 2 . Figure |6] illustrates the 
rapid increase of with /, for selected values of S for 
a rectangular cavity with <r w = 10 6 S/m and fi w<1 — 1. For 
S = 100 m 2 , the predicted values of at / = 0.1, 1 and 
10 GHz are 0.0031, 0.99 and 313, respectively. 

Large values of M may occur, even at relatively low modal 
spectral densities and narrowband excitation, provided the 
spectral overlap of modes is sufficiently high [40] (which 
depends on ohmic dissipation and leakage of the cavity), or 
when the excitation bandwidth is large relative to the average 
spectral modal spacing. For ultra-low loss enclosures (e.g., 
superconducting or laser cavities), the modal overlap can be 
small or nonexistent, even at very high frequencies. Therefore, 
the value of M need not always be exceedingly large, even in 
overmoded conditions. For RF and microwave metal cavities, 

7 Note that Sf is defined as the width at half height. Excitation near the 
edges of the band Sf is not as strong as near the centre of the band but may 
contribute through multiplication of SM C by a suitable factor of order unity. 
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Fig. 6. Moo(/) for a rectangular cavity with er w = 10 6 S/m and 
Mw,r = 1 at selected values of 5*. 



however, the losses and modal overlap are larger whence M is 
then typically large. Since the value of M can also be related 
to the number of coherence cells in V, its value can also be 
estimated by dividing V by the size of such cells [41], [42]. 



